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1.  Introduction 

The  premise  of  this  project  was  that  there  are  non-stationary  signals  that  are  very  difficult  to 
detect  using  conventional  detection  equipment.  This  includes  receivers  which  detected 
correlated  signals,  and  modern  spectrum  analyzers  which  respond  to  voltages  rather  than  the 
peak  signals  of  earlier  detection  signals.  Many  military  listening  devices  assume  that  a  signal 
being  received  is  stationary.  For  example  a  communication  signal  is  a  modulated  carrier  signal 
and  while  not  being  absolutely  stationary  it  can  be  detected  using  a  stationary  assumption  as  on 
average,  that  is  after  a  dwell  interval,  the  signal  can  be  treated  as  stationary.  Other  unintended 
emissions  such  as  a  leaking  clock  signal  are  stationary  and  even  very  low-level  such  signals  can 
be  detected.  Key  to  the  ability  to  examine  such  signals  was  the  development  of  a  scheme  for 
detecting  very  small  signals  close  in  frequency  to  a  large  signal  using  only  power  measurements. 
This  was  based  on  an  earlier  measurement  system  but  enhancements  were  made  to  the 
detection  algorithm  which  better  responded  to  power  than  voltage.  The  system  was  used  to 
measure  a  reflected  radar  signal  from  vibrating  objects  which  has  small  Doppler-like  sidebands. 

This  work  attempted  to  detect  noise  signals  and  chaotic  signals.  We  over  reached,  that  is  for 
sure,  as  detection  of  such  a  signal  proved  particularly  elusive.  We  did  write  two  journal  papers 
from  the  work  performed  under  the  sponsorship  of  this  project  but  much  remains  to  be  done. 
Prior  to  the  project  we  had  established  the  world  record  for  detecting  small  signals  close  in 
frequency  to  a  large  signal.  In  this  project  we  enhanced  this  work  and  in  particular  developed  a 
way  of  measuring  such  signals  even  if  they  were  non-stationary.  The  first  paper  describes  the 
measurement  scheme  and  the  second  paper  applied  the  work  to  the  measurement  of  radar 
reflections  form  a  chaotically  vibrating  antenna.  When  an  antenna  vibrates  it  is  nearly  always  a 
chaotic  vibration.  Or  in  the  preferred  language  of  chaos  investigators  the  antenna  vibrations  have 
non-stationary  nonlinear  dynamic  behavior.  This  description  is  preferred  because  there  are 
many  "chaotic"  behaviors  and  calling  something  chaotic  implies  that  the  behavior  is  reasonably 
well  understood.  The  journal  paper  on  measurements  of  a  chaotically  vibrating  antenna  provides 
a  compact  description. 

This  report  focuses  on  the  results  of  the  research  that  are  not  represented  by  or  not  fully 
described  by  the  papers. 

2.  Current  Signal  Identification  of  Radio  Frequency  Circuits 

This  project  undertook  the  investigation  of  nonstationary  phenomena  at  radio  and  microwave 
frequencies.  The  work  was  based  on  the  premise  that  the  analysis  and  simulation  tools  developed 
for  our  understanding  of  RF  and  microwave  electronics  and  signals  in  general  assume  that  the 
signals  are  stationary  and  that  electronic  phenomena  is  governed  by  integer  calculus.  The  project 
undertook  the  investigation  of  analyzing  signals  and  simulating  electronics  circuits  without 
necessarily  using  integer  calculus.  The  signals  in  circuits  described  by  integer  calculus  are 
deterministic  and  stationary.  Signals  such  as  chaotic  signals,  the  generation  of  which  always 
seems  to  require  nonlinear  feedback  in  a  circuit,  are  deterministic  but  a  not  stationary.  Thermal 
noise  is  a  signal  that  is  neither  deterministic  nor  stationary  but  our  detection  of  noise  assume 
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certain  statistics.  It  is  not  clear  that  we  have  any  ability  to  detect  signals  that  do  not  have  simple 
statistics,  are  non-stationary,  and  are  not  deterministic.  The  research  report  here  made  the  first 
in-roads  but  much  remains  to  be  done. 

Analysis  of  all  virtually  all  electronics,  and  not  just  RF  and  microwave  electronics,  is  based  on  the 
concept  that  signals  are  stationary  and  circuits  act  deterministically.  There  are  a  couple  of 
exceptions.  Very  simple  chaotic  circuits  have  been  considered  previously.  These  are  circuits  that 
generally  have  a  few  reactive  elements  and  just  one  nonlinear  resistive  element.  What  is  known 
as  the  Chua  circuit  is  one  such  circuit.  Most  of  the  investigations  seemed  to  address  how  such  a 
simple  circuit  could  have  such  widely  varying  responses  with  small  changes  in  parameters. 

Our  ability  to  exploit  the  physical  environment  either  through  sensing  or  communications,  for 
example,  is  largely  based  on  partitioning  of  the  physical  space  into  discrete  physical  effects  and 
advances  in  detection  driven  by  signal  processing  enhancements.  Domain  views,  such  as  either 
the  frequency-domain  or  time-domain  views,  are  used  to  simplify  the  effort  required  to 
understand  physical  responses.  There  are  other  views  such  as  the  wavelet-domain  which 
facilitate  the  understanding  of  phenomenology  simultaneously  localized  in  frequency  and  time. 
Domain  views  simplify  a  problem  by  retaining  essential  aspects  and  discarding  nonessential 
aspects.  For  example,  in  the  frequency  domain  view  transient  responses  are  dispensed  with  and 
only  the  steady  state  response  is  retained.  There  may  be  other  domain  views  that  are  ideally 
suited  to  detecting  responses  from  probed  targets  that  could  include  nonlinear  responses  or 
memory-effect  responses. 

Another  example  of  our  inability  to  understand  the  way  the  world  works  is  the  poor  state  of 
speech  recognition.  While  a  speech  recognition  was  not  an  outcome  of  this  research,  it  is  likely 
that  the  way  we  are  going  about  solving  this  particular  problem  is  not  correct.  Could  it  be  that 
we  are  trying  to  force  a  particular  solution  on  a  behavior  which  simply  cannot  be  resolved  using 
our  current  analysis  approaches.  The  brain  surely  does  not  do  arithmetic  with  complex  numbers. 
There  must  be  other  ways  to  process  signals  and  these  should  be  inherently  simple  ways  at  that. 
We  need  different  ways  of  processing  and  extracting  information.  Are  we  making  the  same 
mistake  in  trying  to  force  a  particular  solution  on  electronic  signals  emanating  from  sources  that 
do  not  have  origins  in  well-behaved  integer  calculus  as  well  as  signals  that  are  inherently  non 
stationary. 

3.  Oscillator  Analysis 

The  project  began  with  the  expectation  that  we  could  show  that  phase  noise  on  oscillators  was 
due  to  a  chaotic  process.  The  author  of  this  report  still  believes  that  it  is  but  was  not  able  to 
demonstrate  that  in  this  project.  Our  investigations  were  hampered  by  the  phase  noise  being  so 
low  in  modern  oscillators.  The  circuit  we  were  using  to  provide  experimental  data  was  a  state-of- 
the-art  voltage-controlled  oscillator  centered  at  5.3  GHz.  The  noise  was  too  low  to  capture  using 
a  real-time  oscilloscope,  and  a  belief  was  established  that  frequency-domain  measurements 
perhaps  were  not  valid  in  measuring  the  power  of  a  non-stationary  signal.  The  task  reported  here 
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was  the  development  of  a  relatively  simple  model  of  the  oscillator  that  could  be  used  in  analysis. 
The  5  GHz  Oscillator  is  shown  in  Figure  1.  MATLAB-based  simulations  were  undertaken  and  each 
simulation  required  about  1  week  of  computing  time.  These  simulations  were  performed  as  part 
of  a  ARO  STIR  project  are  reported  elsewhere,  see  the  final  report  for  the  USARO  grant  W911NF- 
15-1-0567. 

This  section  provides  sufficient  information  to  use  the  simplified  model  in  future  research. 

3.1  Extraction  of  Oscillator  Model 

Extraction  of  a  useable  simple  model  of  an  oscillator  from  a  full  blown  circuit  model  necessitates 
an  understanding  of  the  way  the  oscillator  was  designed.  The  majority  of  microwave  oscillators 
are  Colbitts  oscillators  operating  as  admittance  oscillators  with  negative  conductance.  As  the 
oscillation  signal  level  increases  the  net  conductance,  the  sum  of  the  nonlinear  conductance  and 
the  load  conductance  goes  to  zero  as  the  nonlinear  conductance  increases  with  increasing 
amplitude.  Ideally,  to  avoid  multi-frequency  oscillation,  the  active  device  conductance  is 
independent  of  signal  level.  These  characteristics  are  shown  in  Figure  2. 

While  designing  a  stable  voltage  controlled  oscillator,  ideally  we  expect  the  device  parameters  to 
be  dependent  only  on  the  amplitude  of  the  oscillations  as  shown  in  Figure  2  and  the  resonator 
parameters  to  be  dependent  only  on  the  angular  frequency  of  the  oscillations.  The  design  of  a 
reflection  oscillator  involves  in  matching  the  inverse  reflection  coefficient  of  the  active  device  to 
the  reflection  coefficient  of  the  resonator:  l/rd=rr  or  rd*rr=l 

Achieving  this  condition  at  radio  frequencies  over  an  entire  bandwidth  range  is  difficult,  as  at 
higher  frequencies,  the  parasitics  in  the  circuit  are  very  significant  and  complicate 
straightforward  design.  The  oscillator  must  satisfy  the  stability  criterion  over  the  entire  tuning 
range  maintaining  the  output  power  constant  and  phase  noise  low.  This  oscillator  implements  a 
scheme  that  compensates  for  the  effect  of  the  parasitics  by  embedding  a  network  that  modifies 
the  effective  active  device  reflection  coefficient  and  enables  one  port  oscillator  design  techniques 
to  be  used. 
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Figure  1:  C-band  VCO  circuit:  (a)  complete  oscillator  with  the  resonator  network  to  the  left  of  the 
cutaway  line  (x-x)  separated  from  the  active  circuit  to  the  right;  (b)  annotated  resonator  network; 
and  (c)  annotated  active  network.  The  Pi  attenuator  (with  294  Q  resistors  in  the  shunt  legs  and  a 
17.4  Q  series  resistor)  is  between  \/out  and  the  50  Q  bandpass  filter. 
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Figure  2  Variation  of  the  device  parameters  with  respect  to  Frequency  and  Amplitude  in  ideal 
scenario. 

While  designing  negative  resistance  oscillators,  we  use  Kurokawa's  one-port  oscillator  stability 
requirement.  To  apply  this  one  port  criterion,  each  of  the  networks-the  active  device,  resonator, 
and  device  termination  are  characterized  as  one-ports. 

If  the  device  with  admittance  Yd=  Gd+jBd  is  connected  to  a  resonator  of  admittance  Yr=  Gr+jBr, 
the  conditions  required  for  oscillation  can  be  determined  at  equilibrium  as  follows: 

Gd(A)  =  Grico)  Bd(A)  =  Br(ix)) 

where,  A  is  the  voltage  amplitude  and  u  is  the  angular  frequency  at  equilibrium. 


Frequency,  / 


Amplitude  A 


Figure  3:  Intersection  of  the  admittance  loci  of  the  device  and  resonator  network. 

These  conditions  can  be  graphically  represented  by  taking  the  negative  of  the  device's 
admittance  i.e  -Yd(A)=  -Gd(A)-jBd(A)  and  the  locus  of  the  resonator's  admittance  i.e.  Y(ii))= 
Gr(u})+jBr(u})  as  shown  in  Figure  3.  In  order  to  achieve  stable  single  frequency  oscillations,  the 
intersection  of  this  inverse  device  reflection  coefficient  locus  and  the  resonator  admittance  locus 
in  the  complex  plane  must  occur  at  a  single  point.  Stability  issues  will  occur  if  there  are  multiple 
intersections,  which  results  in  multiple  oscillations,  oscillator  start-up  problems  and  unwanted 
noise.  For  an  oscillator  which  operates  over  a  range  of  frequencies,  as  the  circuit  parasitics 
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significantly  affect  the  tuning  behavior,  the  device  parameters  would  be  frequency  dependent. 
Thus,  the  device  admittance  should  now  be  described  as -Yd(A,{j)). 

3.2  Description  of  the  Oscillator 

The  design  objective  of  the  oscillator  in  Figure  4  is  the  generation  of  a  frequency  independent 
negative  conductance  Gd(A,  co)  with  a  prescribed  reflection  coefficient  shape  Td  using  a  transistor 
in  common  base  configuration  with  series  inductive  feedback.  The  effective  device  conductance 
is  modified  when  there  is  reactive  loading,  due  to  which  it  becomes  frequency  dependent. 


Figure  4  Tunable  SiGe  BJT  VCO 

A  negative  conductance  oscillator  is  realized  using  a  series  inductor  in  the  base  lead.  The  base 
inductance  provides  the  required  negative  feedback  between  the  input  and  output  of  the 
transistor  which  results  in  a  negative  conductance  from  the  emitter  to  the  ground.  Since  the 
feedback  is  frequency  dependent,  this  induced  negative  conductance  is  also  frequency 
dependent.  The  reactive  loading  by  the  capacitors  Co  and  Cb,  which  are  chosen  by  reflection 
coefficient  shaping,  modifies  the  effective  device  conductance  so  that  it  now  becomes  frequency 
independent.  The  resonator  network  is  resonant  at  a  frequency  above  the  oscillation  frequency, 
so  it  presents  an  inductance  at  the  oscillation  frequency.  The  resonator  network  consists  of  a 
transmission  line  which  is  coupled  to  stack  of  four  varactors  which  provides  variable  capacitance. 

3.3  Transformation  into  Colpitts  Oscillator 

3.3.1  Theoretical  Explanation 

All  oscillators  can  be  classified  as  either  negative  resistance  oscillators  (reflection  oscillator)  or 
feedback  oscillators.  The  distinguishing  criterion  is  whether  there  is  an  obvious  feedback  circuit 
(a  distinct  feedback  block)  .  At  Radio  frequencies,  it  becomes  difficult  to  construct  a  feedback 
circuit  without  introducing  excess  phase  shift.  Therefore,  almost  all  oscillators  in  the  radio 
frequency  range  are  classified  as  negative  resistance  oscillators. 
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In  a  negative  resistance  oscillator  model,  oscillators  are  represented  as  a  combination  of  an  active 
element  and  a  passive  element.  The  active  element  generates  negative  resistance  and  drives  the 
passive  element,  \A/hich  is  usually  a  resonator.  When  the  negative  resistance  beats  all  losses  in 
the  resonator  as  \A/ell  as  the  active  element,  oscillation  builds  up.  Even  though  the  negative 
resistance  analysis  is  so  common,  still  there  are  questions  regarding  the  validity  of  this  analysis. 
Also,  some  important  oscillation  parameters  like  loaded-Q  cannot  be  derived  from  this  analysis. 
Loaded-Q  is  one  of  the  most  important  fundamental  parameters  because  it  dictates  spectral 
purity  i.e  phase  noise. 

On  the  other  hand,  the  feedback  model  of  an  oscillator  is  known  to  give  better  insights  into 
oscillation  performance.  Although  building  an  obvious  feedback  circuit  is  almost  impossible  at 
radio  frequencies,  interpreting  a  negative  resistance  oscillator  using  the  feedback  model  is 
possible.  This  transformation  enables  transmission  analysis  of  RF  and  microwave  oscillators  that 
are  otherwise  treated  with  the  negative  resistance  model. 


Figure  5:  Replacing  the  resonator  network  with  equivalent  inductance 

In  the  Figure  4  voltage  controlled  oscillator,  we  know  that  the  resonator  network  resonates  at  a 
frequency  above  the  frequency  of  oscillation.  So,  the  resonator  would  present  an  inductance  to 
the  active  network  at  the  frequency  of  oscillation.  Thus,  the  resonator  can  be  replaced  by  an 
equivalence  inductor  as  shown  in  Figure  5  at  the  oscillation  frequency.  Here  the  capacitance  Cx  is 
the  equivalent  capacitance  which  is  obtained  from  the  capacitance  Co  and  base-emitter 
parasitics.  The  Cce  i.e.  the  collector-emitter  capacitance  is  one  of  the  dominant  parasitic 
capacitances  which  contribute  to  the  feedback  network. 


Figure  6  Combining  the  elements. 
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The  inductor  U,  which  is  the  equivalent  of  the  resonator  network,  and  the  equivalent  capacitance 
Cx  can  be  combined  to  form  a  capacitance  Ci  at  the  emitter.  The  collector-emitter  parasitic 
capacitance  can  be  written  as  capacitance  C2.  This  simplification  is  shown  in  the  Figure  6. 

Now,  as  all  the  feedback  elements  Ci,  C2  and  L3  are  obtained,  we  can  replace  the  transistor  by  an 
amplifier  block  and  we  can  represent  this  as  a  feedback  oscillator  as  shown  in  Figure  7.  Thus,  the 
reflection  coefficient  voltage  controlled  oscillator  is  simplified  into  a  basic  feedback  oscillator 
with  amplifier  stage  and  feedback  network. 


Figure  7  Final  feedback  form. 


3.3.2  Derivation  from  transistor  model 

The  above  derived  feedback  form  can  be  extracted  for  the  voltage  controlled  oscillator.  To  start 
with,  we  first  obtain  the  transistor  model  with  all  the  parasitics  listed  in  the  Figure  8. 
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Figure  8:  NE894M13  Si  BJT  transistor  model  with  all  parasitics. 
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From  the  model  we  can  list  the  parasitics  as:  Cce  =  0.05  pF,  Ccb  =  0.01  pF,  Ccbpkg  =  0.05  pF,  Ccepkg  = 
0.05  pF,  Le  =  0.42nH,  Lb  =  0.3  nH,  /.epkg  =  0.05  nH,  /.cpkg  =  0.05nH,  /.bpkg  =  0.05nH.  The  value  of  the 
capacitors  Ca  and  Cb  is  0.5pF  which  are  obtained  from  reflection  coefficient  shaping.  With  the 
VCO  tuned  for  5.3  GHz  oscillation  the  other  parameters  of  the  simple  circuit  are  Ci=  1.8  pF,  C2  = 
0.68  pF,  /.3=2.2nH. 

3.3.3  Summary 

From  the  theoretical  extraction  and  hand  calculations  we  have  obtained  the  values  of  the 
feedback  elements  as  Ci=1.8pF,  C2=0.68pF  and  /.3=2.2nH.  Now  using  these  values  as  starting 
point,  we  can  build  up  the  System  Level  Oscillator  and  then  modify  the  values  accordingly  using 
AWR  simulations  as  yardstick  to  match  this  to  the  Voltage  Controlled  Oscillator. 

3.3.4  Simulation  of  Complete  Oscillator  Model 

The  complete  oscillator  model  agrees  very  well  with  measurements  of  the  fabricated  oscillator. 
In  this  study  the  simplified  model  was  derived  to  match  the  simulated  results  from  the 
performance  of  the  full  oscillator  model.  Both  were  tuned  for  5.3  GHz  operation.  This  subsection 
reports  the  simulated  results  obtained  using  harmonic  balance  simulation  of  he  full  VCO  circuit 
and  the  simplified  model  of  the  oscillator  tuned  for  5.3  GHz  operation.  The  full  oscillator  model 
is  given  in  Figure  9.  Key  simulated  results  are  shown  in  Figures  10  to  13. 
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Figure  9:  Complete  VCO  circuit  schematic  in  AWR. 
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Figure  10:  Tuning  Range  of  the  Oscillator 
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Figure  11:  Spectral  representation  of  the  Oscillator 
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Figure  12:  Output  Current  waveform  for  the  oscillator 
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Figure  13:  Output  Voltage  waveform  for  the  oscillator 
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3.3.5  Oscillator  performance  Using  the  Simplified  Oscillator  Model 

The  simplified  model  in  the  AWR  environment  is  shown  in  Figure  14. 
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Figure  14:  Schematic  of  the  simplified  oscillator  model. 

The  key  simulated  results  using  the  simplified  model  are  shown  in  Figures  16  to  20.  Figure  21 
shows  how  the  quantities  plotted  were  measured. 
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Figure  16:  Spectral  Representation  of  the  simplified  oscillator  model. 
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Figure  17:  Output  Current  waveform  for  the  simplified  oscillator  model. 


Figure  18:  Output  Voltage  waveform  for  the  simplified  oscillator  model. 
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Figure  19:  Transfer  Function  Plot,  Output  Current  vs  Input  Voltage 
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Figure  20:  Determining  the  Frequency  of  oscillation 
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CAP 
ID=C4 
C=0.05  pF 


Figure  21:  Measuring  the  terminal  currents  and  voltages  for  the  transistor 

3.3.6  Summary 

From  the  above  AWR  simulations  of  the  voltage  and  current  \A/aveforms  for  the  Complete 
Oscillator  and  the  System  Level  Oscillator,  we  can  confirm  that  the  system  level  is  a  valid 
extraction  of  the  complete  oscillator.  The  required  delay  in  the  tanh  function  of  the  oscillator  is 
0.125  ns  -  0.0474  ns  =  77.6ps. 
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3.3.7  Oscillator  Simulation  in  MATLAB 

The  simple  oscillator  model  was  simulated  in  MATLAB  with  the  oscillating  waveform  shown  in 
Figure  22. 
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Figure  22:  A  1  GFIz  Colpitt's  oscillator  with  start-up  waveform:  (a)  minimum  circuit;  and  (b) 
simulated  transient  response. 


3.3.8  Phase  Noise  Investigations 

One  of  the  speculations  for  phase  noise  is  that  it  could  be  due  to  time-delayed  feedback  around 
a  weak  nonlinearity  —  this  produces  chaotic  behavior.  A  small  computer  program  was  written 
to  model  the  Colpitt's  osciallator  circuit  shown  in  Figure  23.  The  pi  arrangement  of  capacitors  and 
inductor  ensures  stable  operation  and  the  nonlinearity  of  the  transconductance  is  modeled  so 
that  the  output  current  is  a  tanh  function  of  a  time-delayed  controlling  voltage.  Preliminary 
results  are  shown  in  Figure  24.  These  calculations  were  done  using  quad  precision  (double  double 
precision  in  the  C  programming  language).  Quad  precision  must  be  emulated  and  so  there  is  a 
significant  increase  in  run  time  over  simulation  performed  in  normal  double  precision.  Oscillation 
self-initiates  and  simulations  are  typically  run  for  10  million  time  points  with  the  first  5  million 
discarded  to  remove  transient  effects.  With  so  many  time  points  efforts  are  required  to  minimize 
use  of  memory  and  ensure  that  the  program  remains  in  on-chip  cache.  One  strategy  we  used  was 
avoid  the  use  of  a  fast  Fourier  transform.  Instead  discrete  Fourier  transforms  at  several  discrete 
frequencies  were  performed  by  multiplying  the  calculated  waveform  at  a  particular  time  with  the 
cos  of  the  radian  frequency  times  time  and  summing  the  result.  It  was  necessary  to  perform  the 
Fourier  transformation  in  quad  precision  as  well.  All  this  results  in  a  3  to  5  day  run  time  to  produce 
the  types  of  results  shown  in  Figure  24.  It  takes  considerable  time  to  determine  the  optimum 
parameters  including  the  time  required  to  skip  initial  transients. 
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Figure  24:  Plot  of  the  phase  noise  of  the  Colpitts  oscillator.  Phase  noise  is  normalized  referenced 
to  the  phase  noise  at  100  Hz  offset  from  the  carrier  and  not  strictly  dB  c  (decibels  with  respect  to 
the  carrier). 
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One  of  the  amazing  results  is  the  1// region  of  the  phase  noise  curve  shown  in  Figure  24.  We 
have  tried  other  nonlinear  functions  and  this  slope  is  still  there.  The  phase  noise  level  however 
is  changed  by  varying  the  time  delay.  In  Figure  24  there  is  also  a  hint  that  the  phase  noise  has 
some  higher  slopes  but  many  more  points  are  required  to  do  this.  It  appears  that  there  could  be 
a  l//^  segment  of  the  curve.  Experimentally  phase  noise  has  been  observed  to  have  several  l/ /" 
regions  where  n  ranges  from  1  to  5  but  is  always  an  integer.  Not  all  of  these  regions  are  observed 
however.  The  exciting  thing  will  be  to  see  if  there  are  multiple  l//"  regions  captured  in 
simulation  for  the  same  oscillator.  If  so  then  we  have  a  reasonable  idea  of  the  ultimate  source 
of  phase  noise.  This  could  possibly  lead  to  enhancements  in  our  ability  to  build  low  noise 
oscillators.  If  we  increase  the  number  of  points  by  a  factor  of  10  we  are  looking  at  run  times  of 
20  to  30  days.  Using  a  parallel  computer  does  not  have  any  affect  as  the  program  runs  on  a  single 
core.  We  have  spent  approximately  5  months  investigating  how  the  simulations  can  run  more 
quickly  beginning  with  a  MATLAB  simulation.  Using  double  precision  arithmetic,  which  is 
hardware  based  speeds  analysis  up  by  about  a  factor  of  20  but  there  is  very  little  detail  remaining. 

3.3.9  Summary 

The  corresponding  measured  results  are  shown  in  Figure  25.  The  phase  noise  measured  has  many 
of  the  characteristics  of  the  simple  model  simulated.  Note  that  the  simulated  simplified  model 
does  not  have  a  source  of  noise  but  it  is  known  that  a  weakly  nonlinear  circuit  with  time-delayed 
feedback  is  a  chaotic  oscillator.  The  above  results  seem  surprising  but  we  were  never  able  to 
convince  ourselves  that  the  results  were  nothing  more  than  aliasing.  Indeed  some  of  our 
investigations  showed  that  really  could  be  aliasing.  But  why  are  the  results  so  close.  Could  it  be 
that  the  phase  node  measurement  equipment  (as  well  as  the  MATLAB  simulations)  is  in  fact 
simply  reporting  an  aliasing  result.  This  is  even  more  complicated  by  the  fact  that  the  phase  noise 
could  be  a  chaotic  signal.  Flow  does  a  phase  noise  measurement  apparatus  report  the 
measurement  of  a  chaotic  signal.  All  this  work  was  support  by  almost  1,000  hours  of  simulations 
but  the  more  investigation  done  the  more  confusing  the  whole  situation  became. 
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1kHz  10  kHz  100  kHz  1MHz  10  MHz 

Frequency  offset 

Figure  25:  Phase  noise  measured  at  the  top  and  bottom  of  the  tuning  range  as  well  as  at  5.1  GHz 
where  phase  noise  is  optimum.  Minimum  phase  noise  floor  -116  dBc/Hz  at  1  kHz  offset,  -160 
dBc/Hzat  10  MHz  offset. 
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4.  Investigations  of  Spurious  Signals  introduced  by  a  Vibrating  Antenna 
4.1  Concept 

A  two-tone  signal  traveling  on  a  microstrip  transmission  lines  will  eventually  result  in  third-order 
distortion  although  ostensibly  the  system  is  linear.  What  occurs  is  analogous  to  what  happens 
when  two  tones  applied  to  an  amplifier  as  shown  in  Figure  26.  The  figure  shows  the  spectrum  of 
two  tones  applied  to  the  input  of  an  amplifier  and  the  spectrum  at  the  output  of  the  amplifier 
has  the  original  two  tones  plus  two  additional  smaller  tones  called  the  third-order 
intermodulation  products.  With  PIM  the  third-order  tones  exist  with  what  is  ostensibly  a  linear 
system. 


Figure  26.  An  amplifier  showing  the  spectrum  of  a  two-tone  signal  at  the  input  (on  the  left)  and 
a  spectrum  at  the  output  (on  the  right)  with  two  additional  tones  called  third-order 
intermodulation  produces. 

Earlier  measurements  we  performed  on  transmission  lines  and  passive  terminations  are 
canonically  shown  in  Figure  27.  There  are  three  distinct  regions  A,  B  and  C.  Our  group  was  the 
first  to  be  able  to  report  measurements  corresponding  to  regions  A  and  B.  It  is  now  known  that 
both  regions  are  due  to  electrothermal  effects  with  the  breakpoint  between  regions  A  and  B 
corresponding  to  the  corner  frequency  of  the  lowpass  thermal  system.  Flowever  the  physical 
origins  of  region  C  are  not  known.  We  have  undertaken  measurements  that  seem  to  indicate  the 
region  C  does  not  exist  if  there  is  a  single  metal-on-substrate  transmission  line  system.  We 
speculate  that  it  may  only  occur  when  an  adhesion  metal  is  used.  To  explore  this  we  have  been 
collaborating  with  the  University  of  Surrey  who  have  excellent  metallization  capabilities.  Our 
earlier  attempts  to  create  a  silver  line  on  a  sapphire  substrate  resulted  in  lines  that  were  only 
stable  for  a  few  days  before  delaminating.  This  is  not  surprising  and  this  is  why  adhesion  layers 
such  as  chrome  are  used.  The  adhesion  layer  is  very  thin  but  even  so  we  believe  that  the  different 
work  functions  in  the  two  metal  system  are  enough  to  result  in  distortion  producing  the  region 
identified  as  C  in  Figure  27.  Our  intent  is  to  characterize  the  PIM  performance  of  a  single  metal 
transmission  line  and  compare  it  to  the  characteristics  of  a  transmission  line  with  an  adhesion 
layer.  The  tricky  part  of  this  is  having  a  silver  conductor  on  sapphire  that  does  not  delaminate. 
We  have  been  getting  ready  for  this  experiment  by  increasing  the  stability  and  dynamic  range  of 
our  PIM  measurement  system  and  conducting  experiments  with  trial  structures  produced  by  the 
University  of  Surrey.  Figure  5  shows  a  silver  on  sapphire  test  structure.  We  have  wire-bonded  to 
structures  and  have  done  simple  pull  tests.  After  one  month  no  delamination  was  apparent.  The 
complete  transmission  line  system  is  now  in  fabrication. 
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Figure  27:  Plot  of  passive  intermodulation  distortion  (PIM)  versus  tone  spacing  in  a  two-tone  test. 


4.2.  Papers 

The  work  done  in  this  part  of  the  project  is  documented  in  two  published  papers. 

1.  I.  Kilgore,  J.  Wetherington,  and  M.  B.  Steer,  Algorithmic  strategies  for  high  dynamic  range 
analog  cancellation,"  Government  Microcircuit  Applications  Conf.  (GOMACtech),  March 
2015. 

2.  I.  M.  Kilgore,  S.  A.  Kabiri,  A.  W.  Kane,  M.  B.  Steer,  "The  effect  of  chaotic  vibrations  on 
antenna  characteristics,"  IEEE  Antennas  and  Wireless  Propagation  Letters,  Vol.  15,  2016, 
pp.  1242-1244 

3.  I.  M.  Kilgore,  A.  W.  Kane,  S.  A.  Kabiri,  and  M.  B.  Steer,  Analysis  of  RF  interference  resulting 
from  antenna  vibration,"  Government  Microcircuit  Applications  Conf.  (GOMACtech), 
March  2016. 

The  first  paper  describes  modifications  to  a  small  signal  measurement  capability.  This  describes 
a  measurement  system  for  measuring  small  signals  that  only  requires  knowledge,  i.e. 
measurement,  of  power  and  not  amplitude  and/or  phase.  This  we  believed  was  very  important 
in  measuring  small  signals.  The  second  and  third  papers  describe  application  of  the  measurement 
system  to  signals  reflected  by  an  antenna.  Here  a  900  MHz  signal  was  fed  and  transmitted  by  an 
antenna.  The  small  signal  measurement  system  measured  the  signals  coming  back  from  the 
antenna.  Additional  content  was  added  due  to  a  special  relativity  effect.  Doppler  is  special  case 
of  the  effect.  The  reader  is  referred  to  the  papers  for  a  complete  description.  What  is  shown  in 
the  figures  below  are  photos  of  the  apparatus  used  as  these  photos  are  not  in  the  published 
papers. 
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Figure  28  shows  the  antenna  mounted  on  the  vibration  table.  This  is  a  $150  vibration  table  which 
is  sold  for  high  school  experiments  and  costs  around  $150.  The  table  exhibited  chaotic  vibrations 
as  reported  in  the  cited  paper.  A  laboratory  grade  vibration  table  costing  thousands  of  dollars 
does  not  exhibit  chaotic  vibrations.  The  vibration  table  uses  what  looks  like  a  speaker  coil  which 
was  driven  by  us  at  100  Hz)  to  induce  vibrations  in  the  antenna.  The  antenna  vibrations  were 
measured  by  a  laser  Doppler  vibrometer. 


Figure  28:  Monopole  antenna  attached  to  a  vibration  table. 


The  spectrum  of  the  reflected  RF  signal  from  the  antenna  is  shown  in  Figure  29.  Spurious 
frequencies  are  seen  at  100  Hz  intervals.  Using  a  doppler  model  we  expect  that  these  spurious 
signals  levels  to  roll  off  linearly  with  respect  to  frequency  away  from  the  center  frequency 
however  a  long  tail  was  actually  observed.  There  are  other  indications  that  the  Doppler  effect 
model  is  not  a  complete  model  including  the  asymmetric  response  observed.  Another  possibility 
is  that  additional  harmonics  present  in  vibration  source,  but  we  measured  with  an  LDV  and  the 
vibration  harmonics  are  30  to  40  dB  down. 
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Figure  29:  Signal  "reflected"  from  vibrating  antenna  \A/ith  1  kHz  span.  The  center  frequency  is  at 
915  MHz  and  was  suppressed  through  cancellation.  (Measurement  system  parameters:  Carrier 
frequency  915  MHz;  reference  Level:  -10  dBm;  Signal  generator  po\A/er:  -5  dBm;  OUT  po\A/er:  - 
20.33  dBm;  cancelled  po\A/er:  -62.06  dBm;  Cancellation  achieved:  41.73  dBm;  Measurement  noise 
floor:  -115  dBm;  measurement  system  dynamic  range:  95  dB. 
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Figure  30:  Signal  "reflected"  from  vibrating  antenna  \A/ith  5  kHz  span. 
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Figure  31:  Signal  "reflected"  from  vibrating  antenna  \A/ith  10  kHz  span. 
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Doppler  harmonic  response,  p  =  0  08OCO0 


Figure  32:  Expected  roll  off  using  Doppler  model.  This  used  the  amplitude  of  the  vibrations 
recored  by  the  LDV  at  corresponding  mechanical  vibration  frequencies. 
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Figure  33:  Magnitude  of  reflected  RF  response  at  frequencies  offset  at  harmonics  of  100  FIz. 


Antenna  vibrations  were  also  induced  using  a  studio  speaker  which  provided  more  consistent 
(but  less  powerful)  vibrations  at  lower  frequencies.  See  Figure  34.  The  collection  of  recorded 
responses  are  in  Figures  35  to  39. 


Figure  34:  Acoustically  excited  vibrations. 
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Figure  35:  Measured  "reflected"  spectrum  with  speaker-induced  (acoustic)  vibrations. 
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Figure  36:  Relative  sideband  levels  with  acoustically  excited  vibrations. 
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Figure  37:  Asymmetric  sidebands  of  acoustically  stimulated  vibrations. 
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Single  tone  incident  on  a  vibrating  surface  = 
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Figure  38  :  Frequency  modulation  formulas. 
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Figure  39:  Calculated  response  obtained  assuming  a  single-tone  FM  response.  That  is,  a  single 
tone  incident  on  a  vibratiing  surface. 


5.  Recommendation  for  Follow-on  Work 
5.1  Summary 

Follow-on  work  to  this  project  should  be  less  aggressive,  and  focus  on  developing  investigative 
situations  where  the  nonlinear  dynamic  behavior  is  exaggerated  so  that  it  can  be  captured  by  a 
real  time  oscilloscope  which  of  courses  has  a  reasonably  low  dynamic  range.  The  equivalent 
analog  dynamic  range  of  a  real  time  oscilloscope  is  DR  =  (6.02  x  Q)  (dB).  Where  Q  is  the  number 
of  bits.  For  a  typical  8-bit  real  time  oscilloscope  this  is  48  dB.  This  compares  to  the  dynamic  range 
of  the  frequency  domain  measurement  system  used  with  a  dynamic  range  exceeding  100  dB. 

Much  work  is  left  to  be  done.  There  is  clearly  unusual  behavior,  for  example,  the  paper  on 
vibrating  antennas  could  actually  provide  the  missing  insight  into  way  co-site  interference 
problems  are  so  intractable.  Some  of  the  work  addressed  in  this  report  include  oscillator  phase 
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noise  investigations.  Phase  noise  of  an  oscillator  is  20  to  30  dB  higher  than  can  be  predicted  using 
conventional  analysis.  Our  premise  was  that  perhaps  the  origin  of  phase  noise  was  chaos.  This 
proved  to  be  too  difficult  a  problem  to  look  as  phase  noise  is  small  and  computational  analysis 
required  excessive  analysis  times.  Even  though  this  project  is  over  work  continues  and  here  the 
so-called  Chua  circuit  which  produces  large  levels  of  chaotic  behavior  should  be  used  in  the 
analysis  of  phase  noise.  The  signals  in  this  circuit  can  be  captured  using  a  real  time  oscilloscope 
as  well  as  requiring  modest  simulation  times. 

What  also  remains  to  be  done  is  analysis  of  the  ability  to  detect  non-stationary  signals,  such  as  a 
chaotic  signals,  using  our  current  measurement  equipment  which  uses  a  stationary  assumption. 
One  of  the  issues  that  came  up  in  this  work  is  whether  a  spectrum  analyzer  or  a  phase  noise 
meter  is  capturing  the  power  in  a  nonstationary  signal.  One  line  of  thinking  is  that  is  does  not. 
This  is  a  project  in  its  own  right  and  was  just  part  of  the  project  being  reported  on. 

5.2  Specific  topic  of  future  investigation  are: 

1.  Does  a  modern  spectrum  analyzer  have  aliasing  errors. 

2.  In  spectrum  analyzers  for  20  or  more  years  ago  a  narrow  bandpass  filter,  a  YIG  tuned 
sphere,  was  used  to  pass  only  signals  in  a  very  narrow  band.  A  bandpass  filter  is  a 
correlator  and  so  it  would  treat  noise  differently  than  a  correlated  signal.  What  is  the 
impact  on  the  signal  level  reported. 

3.  Investigate  phase  noise  of  a  noisy  oscillator  with  high-level  noise.  An  example  is  the  Chua 
circuit  that  oscillates  chaotically  and  the  levels  of  chaos  are  large.  Presumably  this  would 
enable  rapid  simulations  and  aliasing  problems  in  the  reported  results  can  be  minimized 
or  at  least  understood.  In  effect  this  is  a  strategy  of  exaggerating  an  effect  so  that  it  can 
be  more  easily  studied. 

4.  Construct  an  RF  noisy  circuit  and  investigate  using  a  high  frequency  real  time  oscilloscope. 

5.  Analyze  the  reflected  signals  from  a  vibrating  antenna  using  a  high-speed  real  time 
oscilloscope. 

6.  Modify  tools  for  the  analysis  of  signals  with  nonlinear  dynamics.  A  chaotic  signal  is  an 
example  of  such  signals.  Current  tools  seem  to  be  focused  on  large-scale  chaotic  effects 
and  proved  not  to  be  robust  for  low  level  signals  in  the  presence  of  large  signals. 
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